Collision system size dependence of dihadron azimuthal correlations in 
ultra-relativistic heavy ion collisions 
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The system size dependence of dihadron azimuthal correlations in ultra-relativistic heavy ion 
collision is simulated by a multi-phase transport model. The structure of correlation functions and 
yields of associated particles show clear participant path-length dependences in collision systems with 
a partonic phase. The splitting parameter (D) and Root Mean Square Width (Afjirms) of away side 
correlation functions increase with collision system size from ^*N-|-^*N to ^^''Au-l-^^'^Au collisions. 
The double-peak structure of away side correlation functions can only be formed in sufficient "large" 
collision systems. These properties provide some hints to study onset of deconfinement, which is 
related to the QCD phase boundary and QCD critical point, by an energy-size scan. 

PACS numbers: 25.75.Gz, 12.38.Mh, 24.85. -|-p 



I. INTRODUCTION 

Quantum Chromodynamics (QCD) calculation pre- 
dicted an exotic quark-gluon matter in QCD phase di- 
agram [l| may be created in the early stage of heavy ion 
collisions at ultra-relativistic energy [31 • For mapping the 
QCD phase diagram and locating QCD phase boundary 
and critical point 0] , one needs to find a way to vary tem- 
perature, T, and chemical potential, /.is. The NA61 col- 
laboration and NA49-future collaboration Q suggested 
that it can be achieved via a systematic energy (E) and 
system size (A) {E — A) scan. The present work focuses 
on the latter in different collision systems with partonic 
phase or hadron gas at ^/snn = 200 GeV. 

Jet quenching phenomenon has been theoretically pre- 
dicted and experimentally observed @. So far, di- 
hadron azimuthal correlation has been demonstrated as 
a good method to reconstruct particle and energy dis- 
tribution induced by the quenched jet. In experiment, 
a double-peak structure was found on the away side of 
dihadron azimuthal correlation functions [9l-[ll| and the 
indication of conical emission of char ged hadrons was re- 
ported by the STAR collaboration A significant 
centrality dependence of double-peak structure of away- 
side correlation functions was observed by the PHENIX 
collaboration and theoretically simulated in Ref . [ll] . 
The centrality dependence results indicate that the struc- 
ture of away-side correlation function is sensitive to sys- 
tem size of the reaction zone. 

These interesting phenomena provoke some theorists 
into explaining the physical mechanisms for origins of the 
double-peak structure, such as Cherenkov-like gluon radi- 
ation mod el 11411 . medium-induced gluon bremsstrahlung 
radiation lla. Il6j. shock wave model in hydrodynamic 
equations |17| , waking the colored plasma and sonic Mach 
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cones [l8|, sonic booms and diffusion wakes in thermal 
gauge-string duality [l^, jet deflection [2^ and strong 
parton cascade mechanism [l3, [2ll - l23 | . 

All of the experimental and theoretical works sug- 
gest that the double-peak phenomenon is a good probe 
to study the hot and dense medium created in ultra- 
relativistic heavy ion collisions 0, HI ■ In this paper we 
investigate the double-peak structure in different collision 
systems, namely ^^N + "N, -f- i^q, ^^-^a + '^^^a, 
27A1 + 27A1, 40ca + 40ca, ^^Cu + 64Cu and i^^Au + 
19'^ Au at = 200 GeV. The collision syste m size 

dependence of the double-peak structure of away-side 
correlation function is focused. We present participant 
path-length i' = 2Nbin/Npart HI (iV(,m and Npart are 
the number of binary collision and participants, respec- 
tively) dependence of the double-peak structure of away- 
side correlation function in the most central collisions 
(0-10%). In addition, the behavior of the away-side cor- 
relation function is also discussed at given various partic- 
ipant path-length v sa 1.90 and 2.28. From these results, 
one can see that the structure of away-side correlation 
function changes near '^'^Ca + ^°Ca collisions at y^SNN 
~ 200 GeV in central collisions (0-10%), but there is no 
significant system dependence of the correlation at given 
participant path-length v. These results show significant 
degree of freedom dependence in the system with a par- 
tonic phase or with a pure hadron gas [26| , which provides 
potential information about the onset of deconfinement. 
This method supports the E—A scan in experiment to in- 
vestigate the QCD phase transition boundary and QCD 
critical point. 

The paper is organized as follows. The following sec- 
tion introduces the model as well as the analysis method 
for dihadron azimuthal correlations. Section III describes 
the detailed results and discussions, which includes the 
structure of dihadron azimuthal correlation functions, 
yields and mean transverse momentum of associated par- 
ticles, the splitting parameter {D) and A(f>r,ns dispersion 
of away-side correlation functions and the corresponding 
results at given participant path-lengths. The nuclear 
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FIG. 1: (Color online) Dihadron azimuthal correlation functions for different collision systems in centrality 0-10% at y^sjviv 
200 GeV; the shadowing area from the STAR data 



modification factor, Rcp, are also presented in Section 
III. Section IV gives a summary. 



II. MODEL AND ANALYSIS METHOD 

In this^ paper, a multi-phase transport model 

(AMPT) |27l, which is a hybrid model, is employed to 
study dihadron azimuthal correlations. It includes four 
main components to describe the physics in relativistic 
heavy ion collisions: 1) the initial conditions from HI- 
JING model [11], 2) partonic interactions modeled by 



a Parton Cascade model (ZPC) [1^, 3) hadronization 
(discussed later) , 4) hadronic rescattering simulated by A 
Relativistic Transport (ART) model [s^ . Excited strings 
from HIJING are melted into partons in the AMPT ver- 
sion with string melting mechanism [3l| (abbr. Hhe Melt 
AMPT version') and a simple quark coalescence model 
is used to combine the partons into hadrons. In the de- 
fault version of AMPT model [3^ (abbr. Hhe Default 
AMPT version')^ minijet partons are recombined with 
their parent strings when they stop interactions and the 
resulting strings are converted to hadrons via the Lund 
string fragmentation model [33|. The Melt AMPT ver- 
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sion undergoes a partonic phase, while a pure hadron gas 
is in the Defauh AMPT version. Details of the AMPT 



model can be found in a review paper j27| and previous 
works [HIIllsl. 

The analysis method for dihadron azimuthal corre- 
lations is similar to that used in previous experiments 
, which describes the azimuthal correlation between 
a high Pt particle (trigger particle) and low pT particles 
(associated particles). The raw signal can be obtained by 
accumulating pairs of trigger and associated particles into 



distributions in the same event. The 



background which is expected mainlyfrom elliptic flow is 
simulated by mixing event method [1, [l^l ■ By selecting 
the associated particles in different events whose centrali- 
ties are every close to that of the same event (raw signals), 
the A(/() distribution can be obtained as the correspond- 
ing background. The background is subtracted from raw 
signal by using A Zero Yield At Minimum (ZYAM) as- 
sumption as that used in experimental analysis [l^ (See 
our detailed analysis in Ref. [T3|). 
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FIG. 2: (Color online) Schematic illustrations, (a) fragments 
from jet in relativistic heavy ion collisions, (b) azimuthal cor- 
relations between trigger particle and low pr fragments. 



III. RESULTS AND DISCUSSIONS 
A. Structure of dihadron correlation function 

The participant path-length, defined as v — 
'^Nhin/Npart |25| , cau describe degree of multiple colli- 
sions between participants in the early stage of heavy 
ion collisions and characterize the size of the reaction 
zone. The n^gp°^^ represents average collision number of 
partons in the Melt AMPT version. The values of and 
significantly increase with varying collision sys- 
tem (CSYS) from "small" size to "large" size at ^/snn = 
200 GeV in the most central collisions (0-10%) as shown 
in Table m From this table, we can see the multiple col- 
lisions are more enhanced in "large" size collision system 
than in "small" size one. The values of Npart, Nhm and 

are comparable to those from the Glauber Model [IHl 
for both 64Cu + 64Cu and ^^^Au + ^^^Au coUisions. 

Figure [1] shows dihadron azimuthal correlation func- 
tions of different collision systems in the most central 
(0-10%) collisions at = 200 GeV. The correla- 

tion functions are calculated in the kinetic windows, 1 
< p^'""" < 3 GeV/c as well as 2.5 < p^"^ < 6 GeV/c 
and jryl < 1. It shows that the structure of away-side 
correlation function changes from the Gaussian-like dis- 
tribution to double-peak structure, obviously emerging 
in '"'Ca -I- 4°Ca collisions, with varying collision system 
from i^'N + i^'N to ^^^Au + ^^^Au coUisions in the Melt 
AMPT version. In this figure, the amplitude of the cor- 
relation function becomes higher with the increasing of 
collision system size. The associated particles in the Melt 
AMPT version are more abundant than those in the De- 
fault AMPT version. 

The yields of associated particles on near-side correla- 
tion functions from the Default AMPT version without 
hadronic rcscattering is very small. The jet pair from HI- 



JING [28[ keeps momentum conservation law and num- 
ber of particles from jet fragmentation dose not balance 
on the near side and away side, which was investigated 
in p + p collisions [IBl- In the Default AMPT version 
without hadronic rcscattering, the Gaussian-shaped frag- 
ments around the jet direction are emitted from the jet, 
while the trigger particle in the fragments are not always 
along the direction of the jet, a schematic illustration 
shown in panel (a) of Fig [5] The azimuthal correlations 
between trigger particle and low pT fragments can be 
obtained in panel (b) in Fig [2] The near-side peak are 
around direction of near-side jet and shifts from direction 
of trigger particle in some cases. As an integral result, the 
amplitude of near side will be very small in the Default 
AMPT version without hadronic rcscattering. 

It is obvious that hadronic rcscattering enhances the 
correlation on the near side and broadens the structure on 
the away side in the Default AMPT version. As we have 
already demonstrated, the structure of dihadron corre- 
lation functions is more reasonable in the Melt AMPT 
version than in the Default AMPT version [31 ■ However 
the Default AMPT version is used to compare the prop- 
erties of the double-peak structure in partonic phase and 
in hadron gas, especially for "small" size collision sys- 
tem. For investigating the properties of collision system 
size dependences of away-side dihadron azimuthal corre- 
lations, we extract the associated particle yield N^^°y, 
splitting parameter (D) (half distance between double 
peaks on the away side ), Root Mean Square Width 
(A</)rms) and mean transverse momentum {pT)away of 
away-side associated particles, which will be discussed 
later, respectively. 
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TABLE I: NpartiCSYS), A^t^^ (CSFS) , u(CSYS) = ^^;^-[g^yf] , in different collision system at ^/i^ = 200 GeV 

for centrality 0-10 %, the value in blanket is taken from Glauber Model [ssj . 
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B. Yield of associated particles 



C. Affirms and splitting parameter on the away side 



The iy{CSYS) dependence of N^^""^ is shown in Fig. [3] 
from the Melt/Default AMPT version, respectively. It 
presents a significant increasing trend of N2^°y with 
varying the collision system from ^""^N + ^^N to ^^^Au 
+ ^^^Au collisions in the most central collisions (0-10%) 
at ./IWn = 200 GeV in the Melt AMPT version. The 
Default AMPT version, with a hadronic gas, does not 
result in a rapidly increasing dependence trend. In the 
Melt AMPT version, the dependence trend indicates the 
jet correlation information can be inherited by more par- 
ticles in a partonic phase than in a hadronic gas, espe- 
cially in "large" size collision system. Furthermore it 
implies that the interaction strength in "large" size colli- 
sion system is more significant than that in "small" size 
collision system, and while strong parton cascade plays a 
dominant role in dihadron azimuthal correlations in the 
Melt AMPT version. It is interesting that the increasing 
slope of N^^°y vs v{CSYS) from the linear fitting in the 
Melt AMPT version is quicker after "^"Ca -f- ''"Ca colli- 
sion system, where clear double-peak structure emerges, 
than that in small systems. This property indicates the 
double-peak (Mach-like) structure can enhance associ- 
ated particles yields of jet correlations partially. 
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FIG. 3: (Color online) Yield of associated particles on away- 
side correlation functions, N^^°y , as a function of v{CSYS) 
by the Melt/Default AMPT version in centrality 0-10% at 
yil^ = 200 GeV. 



Root Mean Square Width {/S.(f)r 
lation function is defined as 
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where A0m is the mean A0 of away-side correlation 
function and it approximates to tt. A(/)„ns can describe 
broadening of away-side correlation function relative to 
TT which is direction of back jet. By A(/)rms, one can in- 
vestigate diffusion degree of the associated particles rel- 
ative to back jet. The v{CSYS) dependences of A^^ms 
in the Melt/Default AMPT version are shown in Fig. H 
respectively. A(f>r,ns from the Melt AMPT version are 
consistent with PHENIX data [13, H^j in Cu Cu and 
Au -f- Au collisions. A^^ms increases from ^^N -1- ^''^N 
coUisions to ^^^Au + ^^'^Au colhsions in the Melt AMPT 
version and the increasing trend is not so quick in the 
Default AMPT version. The increasing trend of A^^ms 
shows broadening of away-side correlation functions with 
increasing size of collision system. It indicates that the 
jet correlation information can reach faraway relative to 
direction of jet with changing the collision system from 
"small" size one to "large" size in a partonic phase. It is 
remarkable that the increasing trend of Acjjrms from the 
linear fitting in the Melt AMPT version shows two differ- 
ent slope after and before "'"Ca + ''"Ca collision system, 
where clear double-peak structure emerges. This inter- 
esting result indicates that the double-peak (Mach-like) 
phenomenon from quenched-jet pushes evolution of jet 
correlations. At given i'{CSYS) case, we investigate sys- 
tem dependence of the double-peak structure in the Melt 
AMPT version. The A(j)rms keeps a flat pattern with 
varying collision system in Fig.[6]at given 1/(0 SYS) « 1.9 
and 2.28. The Acfirms independence of collision system at 
given ^{CSYS) implies that the broadening of the away- 
side correlation function is only sensitive to how violent 
interactions in collision system, which determines how 
far the jet correlation information spreads. These results 
suggest that the back jet modification in the medium 
created in heavy ion collisions with a partonic phase is 
more obvious in the "large" size collision system than in 
"small" size one. 
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FIG. 4: (Color online) Acjtrms of the away-side correlation 
functions as a function of ^{CSYS) in the Melt/Defualt 
AMPT version for the centrality 0-iO% at ^/sj^ = 200 GeV; 
Square from PHENIX data [11 H. 



FIG. 7: (Color online) Mean transverse momentum 
i{PT)a'way) of the away-slde correlation functions as a function 
of u{CSYS) in the Melt/Default AMPT version in centrality 
0-10% at = 200 GeV. 
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FIG. 5: (Color online) Splitting parameter (D) of the away- 
side correlation functions as a function of u{CSYS) in the 
Meh/Defuah AMPT version for the centrality 0-10% at 
./sWn = 200 GeV; Square from PHENIX data 
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FIG. 6: (Color online) A0rms on away-side of the correlation 
functions at given v{CSYS) « 1.90, 2.28 at ^/JWn = 200 
GeV. 



The splitting parameter [D) is another useful observ- 
able to characterize the structure of the double-peak of 
away-side correlation function, and further discloses es- 
sential of jet modification. The v{CSYS) dependence 
of splitting parameter [D) in the Melt/Default AMPT 
version are shown in Fig. [Sj respectively, in the most 
central collisions at ^/snn = 200 GeV. Double Gaussian 
function is not used to fit away-side correlation functions 
since there is no obvious double-peak structure from in 
+ to in 27A1 + 27^1 collisions. The splitting 
parameter (D) increases from "small" size collision sys- 
tem to "large" size one in the Melt AMPT version. It 
is obvious that the splitting parameter {D) is smaller 
in the Default AMPT version than in the Melt AMPT 
version. The results in the Melt AMPT version are com- 
parable to PHENIX data [ll| in Cu -f- Cu and Au + 
Au collisions due to effect of parton cascade in the Melt 
AMPT version [l^ . The parton interaction cross section 
is taken to be 10 mb, which is reasonable for reproducing 
elliptic fiow and dihadron azimuthal correlations in the 
Melt AMPT version [H, [H iH, HI . In addition. Refer- 
ence @ demonstrated that there is a strong suppression 
of inclusive yield and back-to-back correlations at high 
Pt in -'^^^Au 4- ^^^Au collisions relative to that in d + 
^^^Au collisions at ^/sNN = 200 GeV. Theoretical stud- 
ies suggest that the emission angle relative to jet should 
be about 1.23 rad for QGP, 1.11 rad for hadronic gas 
and zero for mixed phase 17[ . The radiation mechanism 
for the double-peak structure [l3] suggests the more en- 
ergetic the jet, the smaller the emission angle. And in 
this calculation it indicates the double-peak structure is 
more obvious in a partonic phase than in a hadron gas. 
A quick change of the double-peak structure is expected 
while one experimentally scan the system size or beam 
energy near the QCD critical point or phase transition 
boundary. It implies that the double-peak structure and 
jet modification is sensitive to the efi'ective degree of free- 
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FIG. 8: (Color online) Rcp in the Default AMPT Version at ^/s]^ = 200 GeV, the shadowing from BRAHMS data [sg] in 
Au + Au collisions. 



dom of the dense medium created in heavy ion colHsions, 
which is related to onset of deconfinement. 

From these results, it can be concluded that a consider- 
able "large" collision system is necessary and the strong 
parton cascade is important for forming the double-peak 
structure of away-side correlation function. An onset of 
the obvious double-peak structure takes place around the 
collision mass range near ^''Ca + ^°Ca collisions. This 
phenomenon indicates that the correlation is sensitive to 
^{CSYS) and n^oi*"", i.e. the correlation depends on 
the collision system size and the violent degree of the 
partonic interaction in a partonic phase. The results in 
different group in a parton phase and a pure hadron gas 
imply the correlations are sensitive to the effective degree 
of freedom of the dense medium created in relativistic 
heavy ion collisions, which can give us some hints of the 
onset of deconfinement in the system-size viewpoint. 



AMPT without parton cascade and hadronic rescatter- 
ing, so we get the high {pT)away this case. And the 
particles will be softened by parton cascade and hadronic 
rescattering. The condition for hadronic rcscattcring is 
sensitive to how close the distance between hadrons is 
in hadronic rescattering (ART) In the Melt AMPT 
version, the system undergoes parton cascade with ex- 
pansion of the system and then partons are combined 
into hadrons. However in the Default AMPT version, 
the system skips the parton cascade stage and most of 
hadrons can participate in hadronic rescattering. There- 
fore hadronic rescattering contribution is more significant 
in the Default AMPT version than in the Melt AMPT 
version. 



E. Suppression of high-pT hadrons 



D. Mean transverse momentum of associated 
particles in central collisions 

For further investigating the contribution of parton 
cascade and hadronic rescattering to jet modification, we 
calculate mean transverse momentum of away-side asso- 
ciated particles in central collisions at y^s^N — 200 GcV 
and is fixed at (0.2, 3) GeV/c. {prhway is pre- 

sented as a function of i'{CSYS) in the Melt/Default 
AMPT version with/ without hadronic rescattering in 
Fig. [3 respectively. {PT)away slightly decreases with 
varying collision system from "small" size to "large" size 
one. {pT)away for different systems are higher in the De- 
fault AMPT version without hadronic rescattering than 
in other case. And {pT)away ^-^c at the bottom and no dif- 
ference in the Melt AMPT/Dcfault AMPT version with 
hadronic rescattering. The results in the Melt AMPT 
version without hadronic rescattering are between them. 
These results are consistent with what we discussed in 
the introduction to AMPT model. The initial condition 
for AMPT is from HIJING [H], corresponding to the 



Jet-quenching can be investigated by dihadron az- 
imuthal correlations as discussed above. Another im- 
portant phenomenon about jet-quenching is suppression 
of high-pT hadrons in relativistic heavy ion collisions. 
The nuclear modification factor, describes this phe- 
nomenon, defined as. 



Rcp 



PTdpTdrj 



[Central) /Ni,i„ (Central) 



d^N 
Pt dpTdrj 



(Peripheral) / Niiin{Peripheral) 



where the central and peripheral collision centralities are, 
respectively, 0-10% and 40-60%. Figure [5] shows the Rcp 
calculated in the Melt AMPT version (right panel) and 
in the Default AMPT version (left panel) in relativis- 
tic heavy ion collisions at ^/snn = 200 GeV. To clearly 
see the system size dependence, we present Rcp in ^^N 
+ i^N, ^^Cu + 64Cu and ^^^Au + ^^^Au collisions and 
those in other collision systems keep the similar depen- 
dent trend. The shadowing area in this figure is from 
BRAHMS data The hadrons in the Melt AMPT 

version, which undergo a partonic phase, are suppressed 
at high Pt and those from the Default AMPT version, in 
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a hadron gas. are not suppressed obviously. The Rcp in 
i^^Au + i^^Au collisions from the Melt AMPT version 
is consistent with the experimental data [11] ■ We can 
conclude that there is jet-quenching-like mechanism in 
the system with partonic phase and jet-quenching phe- 
nomenon is more obvious in the "large" size collision sys- 
tem than in "small" size one. 



cles (-/V^^^^f), A0rms of away-side correlation functions 
and splitting parameter [D) show significant system size 
dependences. These results also present the degree of 
freedom dependence, which is related to onset of decon- 
finement. The present analysis may shed light on the 
QCD phase boundary and critical point in experiment 
by the energy-size scan. 



IV. SUMMARY 

In summary, the present paper discusses the collision 
system size dependences of dihadron azimuthal correla- 
tions at y/SNN = 200 GeV by a multi-phase transport 
model. It presents the properties of away-side correla- 
tion functions, such as N^^J^^, IS.(f)rms, D and {p^""^). 
The onset of obvious double-peak structure occurs near 
"^^Ca -I- ^''Ca collisions and away-side correlation function 
becomes more and more broadening with the increasing 
of collision system size. The yields of associated parti- 
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